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The mechanism responsible for the extraordinary interface conductivity of LaAlOs on SrTiOs and 
its insulator-metal transition remains controversial. Here, using density functional theory calcula¬ 
tions, we establish a comprehensive and coherent picture that the interplay of electronic reconstruc¬ 
tions, lattice distortions, and surface oxygen vacancies fully compensates the polarization potential 
divergence in LaAlOs/SrTiOs, explaining naturally the experimental observations under different 
conditions. While lattice distortions and a charge redistribution between LaO and AIO 2 sub-layers 
play a dominant role in insulating state, a spontaneous appearance of 1/4 oxygen vacancies per 
AIO 2 sub-layer at the LaAlOs surface accompanied by 0.5e“ charge-transfer into the interface is re¬ 
sponsible for interface conductivity and the discontinuous transition in LaAlOs/SrTiOs. Our model 
also explain properties of LaAlOs/SrTiOs prepared with different growth conditions. 


I. INTRODUCTION 

The interface of dissimilar oxide materials hosts rich 
varieties of exotic phenomena not found in its constituent 
materials and has attracted a tremendous amount of re¬ 
search interests, both for fundamental physics and prac¬ 
tical applications^^®. An ideal example is the interface 
of two different insulators, polar LaAlOs on nonpolar 
SrTiOs (LaAlOs/SrTiOs). In their bulk forms, both 
LaAlOs and SrTiOs are wide band gap non-magnetic 
insulators®. Remarkably, when a thin film of LaAlOs 
was deposited on SrTiOs, the interface was found to 
exhibit unusual phenomena such as a conducting two- 
dimensional electron gas (2DEG) with high mobility^’®, 
superconducting® and in some cases irragnetic proper- 
ties^®^^^, and two-dimensional co-existence of both super¬ 
conducting and magnetic properties^®’^^. In particular, 
the fundamental phenomenon that has generated a lot 
of interests is the LaAlOs-thickness-dependent insulator- 
metal transition of LaAlOs/SrTiOs interface. This phe¬ 
nomenon has several peculiar characteristics as shown by 
the various experimental observations summarized in Ta¬ 
ble 1. 

Recent high-energy reflectivity measurements have 
shown that in conducting LaAlOs/SrTiOs, a charge 
transfer of 0.5e“ from LaAlOs into LaAlOs/SrTiOs in¬ 
terface and oxygen vacancies have been observed®. While 
in insulating LaAlOs/SrTiOs, an intra-layer charge redis¬ 
tribution within the LaAlOs film ha-s been found. How¬ 
ever, a fundamental mechanism to explain properties in 
insulating, conducting, as well as an important step- 
function-shape insulator-metal transition as a function 
of LaAlOs thickness in LaAlOs/SrTiOs as one coherent 
picture is still lacking. 

Previous theoretical studies^®^®^ discussed only some 


of experimental observations (in Table I). For in¬ 
stance, calculations based on perfect LaAlOs/SrTiOs 
structures^^’^® or buckling^® did not describe at least two 
very important aspects observed in experiments. First, 
they did not explain the experimentally observed large 
discontinuity in the interface charge density as LaAlOs 
thickness increases above 4 unit cells (uc), but instead 
such calculations predicted this discontinuity only as a 
very gradual change. Second, the classical electronic re¬ 
construction model and previous computational results 
would expect the surface of LaAlOs/SrTiOs to also be 
metallic, which could strongly disrupt any conclusion 
drawn from transport and high-energy optical conduc¬ 
tivity measurements regarding the interface charge den¬ 
sity. Experimental observations instead show that the 
surface is insulating. Furthermore, recent first-principles 
calculations of a polarity-induced defect model®^ were 
not adequate to explain important observations such as 
internal charge redistribution® and lattice distortions®®’^"* 
of LaAlOs in insulating LaAlOs/SrTiOs. 

Here, via first-principles calculations, we propose a 
unified picture that is an interplay of electronic recon¬ 
structions (both between LaAlOs sub-layers, and be¬ 
tween LaAlOs and SrTiOs), lattice distortions, and sur¬ 
face oxygen vacancies compensates the polarization po¬ 
tential divergence introduced due to the polar nature of 
LaAlOs, m different LaAlOs/SrTiOs samples. Using this 
model, we are able to comprehensively explain the funda¬ 
mental mechanism behind the observed physical proper¬ 
ties (see Table I) as one coherent picture. Which of these 
mechanisms are dominant in a given thickness of LaAlOs 
are determined by the lowest total energy state in that 
particular set of conditions. We also show that our model 
can be applied to the LaAlOs/SrTiOs prepared with dif¬ 
ferent growth conditions, and other polar/nonpolar oxide 


2 


TABLE I. A list of fundamental experimental observations 
from various experimental methods of transport, high-energy 
optical conductivity, and structural features. 


Observations for < 4 unit cells (uc) of LaAlOa 

Electrical 

transport 

1. The samples were all insulating. 

High-energy 

optical 

conductivity® 

2. The samples were all insulating. 

3. A redistribution of 0.5e“ from AIO 2 sub¬ 
layer to LaO sub-layer has occurred. 

4. There was no inter-layer charge transfer 
into the interface. 

5. No oxygen vacancies signature has been 
observed. 

Structural 

features^® 

6. An atomic buckling between cations and 
oxygen atoms was observed. 

7. These buckling decreased with increasing 
LaAlOa thickness. 

8. SrTiOa exhibited an opposite buckling 
at the near-interface region, which increased 
with increasing layer thickness. 

Observations for > 4 unit cells (uc) of LaAlOa 

Electrical 

transport®^® 

9. The samples were all conducting. 

10. A sharp discontinuity of the electron gas 
density and insulator-metal transition at 4 
uc LaAlOa has occurred. 

11. If LaAlOa thickness was further in¬ 
creased, the interface conductivity has re¬ 
mained relatively constant, independent of 
LaAlOa thickness. 

12. The surface of LaAlOa was always 
insulating. 

High-energy 

optical 

conductivity® 

13. A decrease of 0.5e“ in LaAlOa him was 
observed. 

14. It was accompanied by an increase of 
0.5e“ at the interface, which resided in the 
SrTiOa side. 

15. Oxygen vacancies signature is observed 
in the LaAlOa him. 

Structural 

features^®’^"* 

16. The buckling for the A1 site for 5 uc 
LaAlOa/SrTiOa him was collapsed. 


electron-ion interaction^^. The cutoff energy for the 
plane wave expansion was set to 400 eV. Gamma cen¬ 
tered k-point grids for Brillouin zone sampling were set to 
3 X 3 X 1 for ionic relaxations and 6 x 6 x 1 for static calcula¬ 
tions, respectively. We used 2x2x1 LaAlOs/SrTiOa su¬ 
percells for both stoichiometric structures and the struc¬ 
tures with oxygen vacancies, the latter of which were cre¬ 
ated by removing one oxygen atom at the surface AIO 2 
sub-layer of the stoichiometric 2x2x1 LaAlOa/SrTiOa 
supercells. The number of SrTiOa layers was fixed to 2 
uc while the thickness of the LaAlOa film increased as 
2, 3, 4, 5, and 6 uc. A much thicker (6uc) SrTiOa as 
substrate tested lead to smilar results reported here. In 
order to minimize the interaction between neighbor sur¬ 
faces, a 13 A vacuum was applied along (001) plane of all 
interface structures. The lateral lattice constant of these 
supercells were fixed to 7.896 A, twice as the equilibrium 
value of the bulk SrTiOa, in good agreement with pre¬ 
vious studies^®’^^. All the atoms except the bottom SrO 
sub-layer were allowed to relax until the forces are smaller 
than 0.02 eV/A. Since SrTiOa is not a polar material, 
the dipole field, if any, must be quite small. To avoid the 
spurious electric field, dipole corrections were included. 
Calculations on the density of states of LaAlOa/SrTiOa 
were performed, confirming that there is no mid-gap state 
being introduced at the bottom of SrO sub-layer. The in¬ 
ternal electric field in LaAlOa was calculated by the slope 
of macroscopically averaged electrostatic potentials^®. 

The formation energy of the surface oxygen vacancies 
is defined by®^: 

Ef = Etotal — Ehost + P)^ 

where Etotal is the calculated total energy of slabs with 
oxygen vacancies, E^ost the energy for the stoichiometric 
2x2x1 supercells, and ijl{T,P) is the chemical potential 
of oxygen, which is related with sample growth environ¬ 
ment. Experimentally, the upper limit of oxygen partial 
pressure for layer-by-layer epitaxial growth of LaAlOa 
films on SrTiOa is 5x10“^ mbar, while the lower limit 
is 1x10“® mbar at the growth temperature of 850 °C^^. 
With these data, we calculate the corresponding energy 
of the upper and lower oxygen chemical potential to be 
-6.77 eV and -7.18 eV, respectively. 


interfaces. 


II. METHODOLOGY 

All calculations were performed by using density- 
functional theory based on Vienna ab initio Simu¬ 
lation Package (VASP)^®’^® with the Perdew-Burke- 
Ernzerhof (PBE) approximation for the exchange- 
correlation functional and the frozen-core all-electron 
projector-augmented wave (PAW) method for the 


III. RESULTS AND DISCUSSIONS 

To investigate the structural and electronic properties 
of stoichiometric 2, 3, 4, 5, and 6 uc LaAlOa layers on 
an SrTiOa substrate, we perform density functional the¬ 
ory (DFT) calculations on the thickness dependence of 
charge density at the interface, averaged lattice distor¬ 
tions of O-La-0 and O-Al-0 chains, and remnant inter¬ 
nal electric field in LaAlOa layers. As shown in Fig. la, 
there is an insulator-metal transition between 4 and 5 uc 
LaAlOs/SrTiOa. For the conducting cases, the charge 
density at the interfaces increases monotonically with 
dLaAiOs, ■ The lattice distortions between O and cationic 
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FIG. 1. Insulator-metal transition and lattice distortions in 
stoichiometric LaAlOa/SrTiOs. (a) LaAlOs thickness depen¬ 
dence of charge density at the interface, along with the aver¬ 
age buckling of La-O-La and O-Al-0 chains in LaAlOs films. 
To study the thickness-dependence behavior, we average the 
angles defined by 180°—a and 180°—/? for O-Al-O and O-La- 
O chains respectively, where a and P are shown in Fig.1(b). 
(b) Structural guide for the buckling of O-Al-O and O-La-O 
chains. We should note that the circled oxygen atoms is not 
in the same plane with La atoms but we define their angle 
P as the projected in-plane angle between their buckling, (c) 
The remnant electric field in LaAlOa films of LaAlOa/SrTiOa 
at different LaAlOa thicknesses. 
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FIG. 2. In-plane averaged charge density of stoichiometric 2 
uc LaAlOs/SrTiOa. (a) Structural configuration of stoichio¬ 
metric 2 uc LaAlOs/SrTiOs and the corresponding charge 
density distribution (isovalue = 0.027 e/A^). (b) Difference 
of the xy-plane-averaged electron density between relaxed and 
unrelaxed cases. 


atoms (as defined in Fig. lb) are found in LaAlOa film 
of all thicknesses. Overall, the lattice distortions of O- 
La-0 chains are more pronounced than the lattice dis¬ 
tortions of O-Al-O chains, consistent with experimental 
observations^^. These cationic displacements introduce 
dipoles opposite to the polar electric field of LaAlOs, and 
thus partially compensates it, leaving a remnant electric 
field in LaAlOg. 

Fig. 1(c) shows the remnant electric field after partial 
compensation by the lattice distortions and electronic re¬ 
constructions. It can be seen that, below 5 uc of LaAlOa, 
the lattice distortions cannot fully compensate the polar 
electric field, as the remnant electric field keeps increasing 
with the increase of d-LaAiOs- When the pLaAiOs > 4 uc, 
a critical point is reached and the remnant electric field 
results in a potential difference that exceeds the LaAlOa 
band gap, which causes an internal Zener breakdown to 
occur. At this point, electronic reconstruction happens in 
LaAlOs/SrTiOs via electrons transfer from surface AIO 2 
sub-layer into interface Ti 02 sub-layers, leaving extra 
holes at surface AIO 2 and extra electrons at interface 
Ti 02 sub-layers, although the calculated charge density 
remains very small. We note that the phenomenon of 
electronic reconstruction was introduced in earlier study 
of polar surfaces in KaCeo^®. This electronic reconstruc¬ 


tion also compensates the polar electric field of LaAlOa, 
leading to a sharp decrease of the remnant electric field 
after 4 uc of LaAlOa (see Fig. Ic). 

Recent high-energy optical conductivity measurement 
of LaAlOa/SrTiOa shows that there is charge redistri¬ 
bution between LaO and AIO 2 sub-layers in insulating 
LaAlOa/SrTiOa®. In Fig. 2, we compare the charge den¬ 
sity distribution between relaxed and unrelaxed stoichio¬ 
metric 2 uc LaAlOa/SrTiOa structures and the internal 
charge redistribution between LaAlOa sub-layers is in¬ 
deed found in our calculations. Because of the lattice 
distortions, the La atoms move upwards to the surface 
while the oxygen atoms in AIO 2 sub-layers move down¬ 
wards to the interface (the La-0 buckling dipole), mak¬ 
ing their orbitals more overlapped and thus more cova¬ 
lent than in the structures without lattice distortions. 
Although the overlap between A1 atoms in AIO 2 sub¬ 
layers and the oxygen atoms in LaO sub-layers (the Al- 
O buckling dipole) counteracts the charge redistribution 
caused by the La-0 buckling dipole, this effect is inferior 
due to the less pronounced lattice distortions of 0-Al- 
O chains than the O-La-O chains. Since this internal 
LaAlOa charge redistribution has strong interplay with 
the lattice distortion effects, it also disappears once the 
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FIG. 3. Formation energy of surface oxygen vacancies in 
LaAlOa/SrTiOa with different oxygen chemical potential and 
LaAlOa thicknesses, (a) Oxygen chemical potential depen¬ 
dence of formation energy of surface oxygen vacancies in 
LaAlOa/SrTiOs. The vertical black dash line indicates the 
oxygen rich condition, while the vertical red dash line indi¬ 
cates the oxygen poor condition, (b) The formation energy 
of surface oxygen vacancies for experimental oxygen-rich and 
oxygen-poor conditions as a function of LaAlOs thicknesses 
dependence in LaAlOs/SrTiOs. 


atomic lattice distortions vanish. For this stoichiometric 
2 uc LaAlOa/SrTiOa structure, there is no charge redis¬ 
tribution between LaAlOa and SrTiOa, consistent with 
experimental observations®. 

Up to this point, we are able to explain to some 
extent the experimental observations of insulating 
LaAlOa/SrTiOa (2 and 3 uc LaAlOa/SrTiOs), i.e. the 
electronic redistribution inner LaAlOg sub-layers®, the 
lattice distortions^®’^^, and the remnant electric field®®. 
However, this stoichiometric LaAlOa/SrTiOa model is 
not able to fully explain the conducting LaAlOa/SrTiOa. 
First, the critical LaAlOa thickness happens at 5 uc in¬ 
stead of 4 uc as observed experimentally®. Second, in 
conducting samples, the charge density increases with 
LaAlOs thickness, while experimentally the charge den¬ 
sity of conducting LaAlOs/SrTiOs is independent from 
LaAlOs thickness, as observed in the step-like insulator- 
metal transition®’®. Third, the charge density is much 
less than the 0.5e“ observed experimentally®. Fourth, 
the lattice distortions still remain after the interface be¬ 
comes conducting, while experimentally^®’^^ they col¬ 
lapse in conducting LaAlOs/SrTiOs. 

Another important observation in high-energy optical 
conductivity of LaAlOs/SrTiOs measurement was that 
the conducting LaAlOs/SrTiOs exhibited Vq signature 
in LaAlOs.® However, this signature was not observed in 
insulating LaAlOs/SrTiOs. This suggests that the for¬ 
mation conditions for Uo in insulating and conducting 
LaAlOs/SrTiOs are different. For this reason, we cal¬ 
culate formation energy (Ey) of Vq in LaAlOs film of 
LaAlOs/SrTiOs with different LaAlOs thicknesses. First 
of all, our calculations show that in LaAlOs film, the E/ 
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FIG. 4. Insulator-metal transition and lattice distortions in 
LaAlOs/SrTiOs, after the introduction of surface oxygen va¬ 
cancies in conducting LaAlOs/SrTiOs. (a) Structural guide 
for the relaxed 4 uc LaAlOs/SrTiOs along < 010 > direction 
with surface oxygen vacancy, (b) Structural guide for the 
relaxed 4 uc LaAlOs/SrTiOs along < 100 > direction with 
surface oxygen vacancy, (c) LaAlOs thickness dependence 
of number of excess charge at the interface, and the average 
buckling of La-O-La and O-Al-O chains in LaAlOs thin films 
for different thicknesses of LaAlOs on SrTiOs. The 2 and 3 
uc LaAlOs/SrTiOs are stoichiometric, while the 4, 5, and 6 
uc LaAlOs/SrTiOs are with surface oxygen vacancies. 


of Vo is the lowest at the surface . Secondly, the E/ 
of surface Vo in LaAlOa/SrTiOa also varies as a func¬ 
tion of LaAlOs thickness and oxygen chemical potential, 
as shown in Fig. 3. One clearly sees that E/ of sur¬ 
face Vo decreases as the thickness of LaAlOs increases, 
which means it is more energetically favorable to cre¬ 
ate surface Vo in thicker LaAlOs film. Interestingly, 
when dLaAlOs, — ^ uc, the formation energy of surface 
Vo crosses through zero at most oxygen chemical po¬ 
tential window, indicating a high possibility to form the 
surface Vo. 

Experimentally, the upper limit of oxygen partial 
pressure for layer-by-layer deposition of LaAlOs/SrTiOs 
is 5x10“^ mbar with a typical deposition temperature 
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of 850 °C^^, implying this oxygen partial pressure is the 
oxygen-rich limit of LaAlOa/SrTiOa deposition. The cal¬ 
culated formation energy of surface oxygen at this oxygen 
chemical potential energy is shown in Fig. 3b. It can be 
seen that formation energy is equal to or below zero for 
d-LaAiOs — 4 uc, even at this oxygen-rich limit, suggest¬ 
ing that the oxygen vacancies are energetically favorable 
to be formed in LaAlOa film with these thicknesses. 

Based on this reason, we perform DFT calculations by 
including surface Yq for 4, 5, and 6 uc LaAlOa/SrTiOa, 
which are experimentally found to be conducting®. The 
relaxed LaAlOa/SrTiOa structure after the inclusion of 
surface Vo is shown in Fig. 4a, b along different direc¬ 
tions. From here, the dipole-inducing lattice distortions 
as defined in Fig. lb vanish, consistent with experimen¬ 
tal results^®’^^. The lattice distortions seen in Fig. 4a 
and 4b are mainly caused by the absence of one surface 
oxygen atom in 2x2x1 supercell, and do not give rise to 
cationic dipoles. 

The results of the calculations for the charge den¬ 
sity are also shown in Fig. 4c. (Note that the calcula¬ 
tions for the charge density and buckling in 2 and 3 uc 
LaAlOa/SrTiOa are performed without including surface 
Vo, the same as the results in Fig. la.) It can be seen 
that the LaAlOa thickness dependence of the interface 
charge density now forms a step function, resulting in a 
discontinuity transition, consistent with experimental ob¬ 
servations®’®. There are four main effects of surface Vo- 
First, the charge density is found to be 0.5e“, which com¬ 
pletely compensates the polar electric field, consistent 
with the experimental observation®. Second, the amount 
of 0.5e“ is consistent with the inter-layer charge transfer 
stipulated by the electronic reconstruction model. Third, 
the 0.5e“ charge density is independent of LaAlOa thick¬ 
ness (and results in an insulating surface), also in good 
agreement with experimental results®’®. Fourth, we can 
reconcile the presence of surface Vo with the reconstruc¬ 
tion model for polar crystal terminations, as discussed 
below. 

In the conventional electronic reconstruction model^®, 
without allowing for Vo or the internal lattice distortions 
described above, the charge of 0.5e“ at the interface is 
compensated by 0.5 holes in the mainly 0-2p band of 
the A102-terminated surface of LaAlOa. In this model 
not allowing for lattice distortions or other deviations the 
electronic reconstruction would occur immediately after 
even only one layer of LaAlOa was deposited (see Fig. 5a), 
contrary with experimental observations®’®. Even more 
importantly, both the SrTiOs interface layer and the top 
AIO 2 surface layer would become metallic. Allowing for 
the lattice distortions, as shown in Fig. 5b, demonstrates 
that these internal distortions can compensate for the 
internal polarization potential divergence for dLaAiOs < 
4 uc, without electronic reconstruction into the interface. 
Meanwhile, for dLaAio^ > 4 uc the top of the valence 
band of LaAlOa and the bottom of the conduction band 
of SrTiOa starts to overlap. These result in both the 
interface and the surface becoming metallic with a very 


a 



b 



c 



5 Number of LAO Layers on STO (uc) 

FIG. 5. A discontinuity phase-like transition due to interplay 
of electronic reconstructions, lattice distortions, and surface 
oxygen vacancies. Three cases of LaAlOs/SrTiOa are consid¬ 
ered: (1) LaAlOs/SrTiOa without both lattice distortions and 
surface oxygen vacancies (blue squares), (2) LaAlOs/SrTiOs 
with lattice distortion but without surface oxygen vacancies 
(red circles), and (3) LaAlOs/SrTiOs with both lattice dis¬ 
tortion and surface oxygen vacancies (black triangles), (a) 
Comparison of the interface charge density of each case, (b) 
Comparison of O-Al-0 chains buckling in LaAlOs film of each 
case, (c) Comparison of relative total energy of the second 
and third cases, with the first case as the reference. Below 4 
uc LaAlOs film (red-shaded region), the second case results 
in the lowest total energy and is consistent with experimental 
observations®’®’^®. On the other hand, above 4 uc LaAlOs 
film (green-shaded region) it is the third case that results in 
the lowest total energy and agreement with experimental ob¬ 
servations®’®’^®. 
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small but gradually increasing charge density of opposite 
sign as the LaAlOa thickness increases (see Fig. la and 
5a). 

The energy to create Vo, however, approaches zero at 
4 uc of LaAlOa (Fig. 3) and becomes negative as LaAlOs 
thickness is further increased. In this case, the total en¬ 
ergy of the system with 1/4 Vo per AIO 2 sub-layer and 
reduced lattice distortions is lower than the total energy 
for the case without the surface Vo and with the strong 
lattice distortions. Thus, this becomes the energetically 
favorable situation (see Fig. 5c). This results in 0.5e“ at 
the interface, while the top AIO 2 sub-layer remains insu¬ 
lating because the electrons left behind by the formation 
of Vo are transferred to the interface. In this way the 
sample as a whole remains charge neutral as must be the 
case. 

Based on these results, our model can naturally explain 
the LaAlOs-thickness-dependent discontinuous transi¬ 
tion at 4 uc LaAlOs of LaAlOa/SrTiOa. Below 4 uc 
of LaAlOa, the internal polarization potential in LaAlOa 
is compensated by the lattice distortions of LaAlOa film 
without electronic reconstruction into the interface ac¬ 
companied by the intra-layer charge redistribution within 
the LaAlOs film, resulting in insulating LaAlOs/SrTiOa. 
Above 4 uc of LaAlOa, the internal polarization poten¬ 
tial in LaAlOs is instead compensated by electronic re¬ 
construction into the interface, stabilized by the ener¬ 
getically favorable spontaneous appearance of Vo in the 
surface AIO 2 sub-layer of LaAlOa, resulting in conduct¬ 
ing LaAlOa/SrTiOa interface with LaAlOa-thickness- 
independent charge density. 

Furthermore, we can classify three types of 
LaAlOs/SrTiOs samples based on different experi¬ 
mental conditions®’^^^^^. First, when the LaAlOa film is 
deposited on SrTiOs with moderate oxygen partial pres¬ 
sure, surface Vq forms, which leads to the insulator to 


metal behavior of LaAlOa/SrTiOa described in Fig. 4c. 
Second, when this LaAlOa/SrTiOa is post-annealed 
in high oxygen pressure, the experimentally-observed 
conductivity decreases, but is not completely gone®’®®’®'^. 
Based on our model, this suggests that the oxygen 
annealing removes some of the surface Vq, leading to 
the insulator-metal transition behavior of stoichiometric 
LaAlOa/SrTiOa described in Fig. la with its much 
lower interface charge density. Third, when LaAlOa 
is deposited with very low oxygen pressure, it is ex¬ 
perimentally observed that the conduction region of 
LaAlOa/SrTiOa extends up to 500 /rm®^. This suggests 
that the low oxygen pressure condition generates a high 
number of Vq in SrTiOa (instead of surface LaAlOs), 
giving rise to three-dimensional electron gas at the 
interface®^. 

In conclusion, our result shows that interplay of elec¬ 
tronic reconstructions, lattice distortions, and surface 
oxygen vacancies is responsible for the insulator-metal 
transition of LaAlOa/SrTiOa with a step function at four 
LaAlOa unit cells and can be controlled by the sample 
growth conditions and environments. 

ACKNOWLEDGMENTS 

We acknowledge Ilya Elfimov, Alex Zunger, and 
Matthew Sherburne for useful discussions and com¬ 
ments on teh manuscript. This work is supported 
by Singapore National Research Foundation under 
its Competitive Research Funding (NRF-CRP 8-2011- 
06 and NRF2008NRF-CRP002024), MOE-AcRF Tier-2 
(MOE2010-T2-2-121), NUS-YIA, and FRC. We acknowl¬ 
edge the NUS Graphene Center and CSE-NUS comput¬ 
ing centers for providing facilities for our numerical cal¬ 
culations. The work at UBC is supported by the Cana¬ 
dian funding agencies NSERC, and CRC as well as the 
University of British Columbia. 


* phyandri@nus.edu.sg 

^ H. Y. Hwang, Y. Iwasa, M. Kawasaki, B. Keimer, N. Na- 
gaosa, and Y. Tokura, Nat. Mater. 11 , 103 (2012). 

^ C. Cen, S. Thiel, G. Hammerl, C. W. Schneider, K. E. 
Andersen, C. S. Hellberg, J. Mannhart , Nat. Mater. 7, 
298 (2008). 

® C. Cen, S. Thiel, J. Mannhart, and J. Levy, Science 323, 
1026 (2009). 

^ P. Irvin, Y. Ma, D. F. Bogorin, C. Cen, C. W. Bark, C. 
M. Folkman, C.-B. Eom, and J. Levy, Nat. Photon. 4, 849 
( 2010 ). 

® N. Nakagawa, H. Y. Hwang, and D. A. Muller, Nat. Mater. 
5, 204 (2006). 

® T. C. Asmara, A. Annadi, I. Santoso, P. K. Gogoi, A. 
Kotlov, H. M. Omer, M. R. Motapothula, M. B. H. Breese, 
M. Rbhausen, T. Venkatesan, Ariando, A. Rusydi, Nat. 
Commun. 5, 3663 (2014). 

^ A. Ohtomo and H. Y. Hwang, Nature 427, 423 (2004). 

® S. Thiel, G. Hammerl, A. Schmehl, C. W. Schneider, and 
J. Mannhart, Science 313, 1942 (2006). 


® N. Reyren, S. Thiel, A. D. Caviglia, L. F. Kourkoutis, 
G. Hammerl, C. Richter, C. W. Schneider, T. Kopp, A.- 
S. Riietschi, D. Jaccard, M. Gabay, D. A. Muller, J.-M. 
Triscone, J. Mannhart, Science 317, 1196 (2007). 

A. Brinkman, M. Huijben, M. Van Zalk, J. Huijben, U. 
Zeitler, J. C. Maan, W. G. Van Der Wiel, G. Rijnders, 
D. H. A. Blank and H. Hilgenkamp, Nat. Mater. 6, 493 
(2007). 

Ariando, X. Wang, G. Baskaran, Z. Q. Liu, J. Huijben, J. 

B. Yi, A. Annadi, A. Roy Barman, A. Rusydi, S. Dhar, 
Y. P. Feng, J. Ding, H. Hilgenkamp, T. Venkatesan, Nat. 
Commun. 2, 188 (2011). 

J.-S. Lee, Y. W. Xie, H. K. Sato, C. Bell, Y. Hikita, H. Y. 
Hwang and C.-C. Kao, Nat. Mater. 12, 703 (2013). 

L. Li, C. Richter, J. Mannhart, and R. C. Ashoori, Nat. 
Phys. 7, 762 (2011). 

J. A. Bert, B. Kalisky, C. Bell, M. Kim, Y. Hikital, H. Y. 
Hwang and K. A. Moler, , Nat. Phys. 7, 767 (2011). 

R. Pentcheva and W. E. Pickett, Phys. Rev. Lett. 102 , 
107602 (2009). 



7 


M. Stengel, and D. Vanderbilt, Phys. Rev. B 80, 241103(R) 
(2009). 

W.-J. Son, E. Cho, B. Lee, J. Lee, and S. Han, Phys. Rev. 
B 79 , 245411 (2009). 

C. Cancellieri, D. Fontaine, S. Gariglio, N. Reyren, A. D. 
Caviglia, A. Fete, S. J. Leake, S. A. Pauli, P. R. Willmott, 

M. Stengel, P. Ghosez, J.-M. Triscone, Phys. Rev. Lett. 
107 , 056102 (2011). 

L. Zhang, X.-F. Zhou, H.-T. Wang, J.-J. Xu, J. Li, E. G. 
Wang, and S.-H. Wei, Phys. Rev. B 82, 125412 (2010). 

N. C. Bristowe, P. B. Littlewood, and E. Artacho, Phys. 
Rev. B 83, 205405 (2011). 

Y. Li, S. N. Phattalung, S. Limpijumnong, J. Kim, and J. 
Yu, Phys. Rev. B 84, 245307 (2011). 

L. Yu and A. Zunger, Nat. Gommun. 5, 5118, (2014). 

S. A. Pauli, S. J. Leake, B. Delley, M. Bjorck, G.W. 
Schneider, C. M. Schlepiitz, D. Martoccia, S. Paetel, J. 
Mannhart, and P. R. Willmott, Phys. Rev. Lett. 106 , 
036101 (2011). 

M. Salluzzo, S. Gariglio, X. Torrelles, Z. Ristic, R. D. Ca¬ 
pua, J. Drnec, M. M. Sala, G. Ghiringhelli, R. Felici, and 

N. B. Brookes, Adv. Mater. 25, 2333 (2013). 


G. Kresse and J. Hafner, Phys. Rev. B 47 , 558(R) (1993). 
G. Kresse and J. Hafner, Phys. Rev. B 49 , 14251 (1994). 
G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 
A. Balderesch, S. Baroni, and R. Resta, Phys. Rev. Lett. 
61 , 734 (1988). 

R. Hesper, L. H. Tjeng, A. Heeres, and G. A. Sawatzky, 
Phys. Rev. B 62 , 16046 (2000). 

M. Takizawa, S. Tsuda, T. Susaki, H. Y. Hwang, A. Fuji¬ 
mori, Phys. Rev. B 84 , 245124 (2011). 

J. Osorio-Guillen, S. Lany, S. V. Barabash, and A. Zunger, 
Phys. Rev. Lett. 96 , 107203 (2006). 

M. Basletic, J.-L. Maurice, C. Carretero, G. Herranz, O. 
Copie, M. Bibes, E. Jacquet, K. Bouzehouane, S. Fusil, A. 
Barthelemy, Nat. Mater. 7 , 621 (2008). 

C. Cancellieri, N. Reyren, S. Gariglio, A. D. Caviglia, A. 
Fete, J.-M. Triscone, EPL 91,17004 (2010). 

Z. Q. Liu, C. J. Li, W. M. Lii, X. H. Huang, Z. Huang, S. 
W. Zeng, X. P. Qiu, L. S. Huang, A. Annadi, J. S. Chen, 
J. M. D. Coey, T. Venkatesan, Ariando, Phys. Rev. X. 3, 
021010 (2013). 


